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C^; ABSTRACT 
O ■ 

^ . The non-Gaussian contribution to the intrinsic halo spin ahgnments is ana- 

'— lytically modeled and numerically detected. Assuming that the growth of non- 
CSI . Gaussianity in the density fluctuations caused the tidal fleld to have nonlinear- 

^ . order effect on the orientations of the halo angular momentum, we model the 

Q\ . intrinsic halo spin alignments as a linear scaling of the density correlations on 

^ . large scales, which is different from the previous quadratic-scaling model based 

[>-^ . on the linear tidal torque theory. Then, we analyze the halo catalogs from 

^ . the recent high-resolution Millennium Run simulation at four different redshifts 

O ■ (2; = 0, 0.5, 1 and 2) and measure quantitatively the degree of the nonlinear effect 

>■ , on the halo spin alignments and its changes with redshifts. A clear signal of spin 

^ . correlations is found on scales as large as 10/z~^Mpc at z = 0, which marks a 

detection of the nonlinear tidal effect on the intrinsic halo alignments. We also 
investigate how the nonlinear effect depends on the intrinsic properties of the 
halos. It is found that the degree of the nonlinear tidal effect increases as the 
halo mass scale decreases, the halo speciflc angular momentum increases, and 
the halo peculiar velocity decreases. We discuss implication of our result on the 
weak gravitational lensing. 



Subject headings: cosmology: theory — large-scale structure of universe 
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INTRODUCTION 



The intrinsic galaxy alignment refers to a cosmic phenomenon that the orientations of 
the galaxy spin axes (or galaxy ellipticities) are not random but locally aligned between 
neighbors due to the spatial correlations of the initial tidal field. The topic of the intrinsic 
galaxy alignments has recently drawn considerable attention in the field of the large-scale 
structure for three reasons. Firstly, since it is believed to be induced by the spatial correla- 
tions of the initial tidal field, it will be possible in principle by measu ring the galaxy in t rinsic 



alignments to reconstruct t he initial density field on galactic scales ( iLee fc Perul200d . 12001 
2002l : lLee fc Erdogdul[2007[ ). 



Secondly, the initial tida l correlations are also responsible for the web-like distribution 
of the large scale structures (iBond. Kofman. fc Pogosyanlll996( ). Thus, the galaxy intrinsic 
alignments will be useful to characterize and quantify the influence of the large-scale struc- 



tures on the orientations of the galaxies distributed in a cosmic web (iNavarro et. al.l 12004 



Truiillo et. al.ll2006l : IPatiri et al.ll2006l : lAraeon-Calvo et al.l 120071 : iHahn et al.l 120071 ) 



Thirdly and most importantly, the galaxy intrinsic alignments could cause non-negligible 
contamination on weak lensing signals. Plenty of efforts have been made so far to distinguish 
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For these reasons, it is of fundamental importance to predict the degree of galaxy in- 
trinsic a lignments with a valid theory and measure quantitatively their signals in practice. 
In fact, iPen et al.l (120001 ) have provided an analytic mo del for the intrinsic galaxy spin 



alignments in the frame of the linear tidal torque theory (IDoroshkevichl Il970l ; IWhitd Il984 



Catelan fc Theunslll996l ). According to their model, the galaxy intrinsic alignments can be 
expressed in terms of a quadratic scaling of the density correlation function and thus are 
expected to exist only at small distances of order of a few Mpc. A crucial implication of their 
analytic model is that due to the quadratic scaling nature of the galaxy intrinsic alignments 
the cross-correlations between the gravitational lensing shears from large scale density field 
and the galaxy intrinsic alignments with local tidal shears will be zero. 



Pen et al.l (I2OOOI ) tested their analytic predictions against low-redshift samples of spiral 
galaxies and found that the observed signals are consistent with the analytic predictions. Yet, 
their detections suffered from poor-number statistics and thus were still tentative. Later, 
several other authors have confirmed the existence of galaxy intrinsic align ments at small 
scales either in numerical simulations or in low-redshift observational surveys (IHeavens et al. 
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2000l : Ijindbooa ISrown et al.ll2002l : iHevmans et alJl2004f ) 



However, iHui &: Zhana (120021 ) has pointed out that it might be invahd to describe the 
galaxy intrinsic ahgnments on large scales as a quadratic function of the density correlation. 
Their logic is as follows: Since the quadratic scaling is based on the linear tidal torque the- 
ory which adopts a somewhat oversimplified assumption that the tidal field is Gaussian in 
the subsequent evolutionary stages, it should not be a good approximation to describe the 
galaxy intrinsic alignments on large scales. In reality, the density fluctuations will develop 
non-Gaussianity via gravity which would in turn lead to non-negligible contributions of the 
nonlinear-order of the tidal tensors to the generation of the galaxy angular momentum. Due 
to this nonlinear effect on the galaxy angular momentum, their intrinsic spin alignments 
should be better approximated as linear scaling with the density correlation function. Be- 
cause the linear scaling of the density correlation drops much slowly than the quadratic 
scaling, the intrinsic spin alignments would not be completely negligible even on large scales. 
If their claims turn out to be true, then it will have a signiflcant impact not only on the 
weak lensing analysis but also on our fundamental understanding of the evolution of the 
tidal alignments. 

In the light of their claims, the following questions naturally arise; Does the nonlinear 
tidal effect on the galaxy intrinsic alignments really exist to a non negligible degree? If so, 
at what epochs and on which scales does its contribution begin to be signiflcant? Does it 
depend on the intrinsic properties of the galaxies or dark halos? Our goal here is to answer 
the above questions using both analytical and numerical methods. 

The outline of this paper is as follows: §2, we overview briefly the previous analytic model 
for the galaxy spin-spin alignments based on the linear tidal torque theory and propose a 
new model to account for the nonlinear tidal effect. §3, we report a detection of the signals 
of the nonlinear tidal effect on the intrinsic spin alignments of dark matter halos simulated 
in N-body experiments and show how the signals depend on redshift, scale, and the halo 
intrinsic property. §4, we summarize our results and discuss the implications of our work on 
the weak lensing effect. 



ANALYTIC MODEL 



2.1. Overview of the Linear Tidal Torque Theory 



The linear tidal torque theory explains that unless a proto-halo has a perfectly spherical 
shape, it can acquire spin angular momentuni at flrst order t hrough its tidal interactio n with 
the surrounding matter (IDoroshkevichl Il970l : IWhitd 1 19841 : ICatelan fc Theund Il996l ) . The 
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main prediction of the linear tidal torque theory is that the proto-halo angular momentum 
vector, (Jj), is proportional to the anti-symmetric product of the two tensors, the inertial 
momentum tensor (lij) and the local tidal shear tensor, (T^), as 



Here the defintions of 7^ and Tij are given as 

lij = / d^qp{q)qiqj, Ti 
Jv 



dqidq. 



-$(q) 



(2) 



where q is the Lagrangian position of the particles that reside in the proto-halo regions, V 
and p(q) are the Lagrangian volume and the density of the proto-halo region, respectively, 
and $(q) is the velocity perturbation potential. 

Given the property of the perfectly anti-symmetric tensor eijk in equation ([1]), an ad- 
ditional condition has to be satisfied for the first-order generation of the proto-halo angu- 
l ar momentum: the prin cipal axes of lij and Tij have to be misaligned with each other 
( ICatelan &: Theunslll996l ). Furthermore, a crucial implication of equation ([T]) is that if this 
additional condition is satisfied and thus the angular momentum of a proto-halo is generated 
at first order, then the direction of the proto-halo angular momentu m is not randorn but 
preferentially aligned with the principal axes of the local tidal tensor (ILee fc PenI 120001 ) . 



Numerical experiments have re vealed that the principal axes of ) and (Tij) are cor- 
related strongly but not perfectly (ILee fc Peru l2000l : iPorciani et al.ll2002l ). which indicates 
that the first-order generation of the proto-halo angular momentum is not so efficient and 
the degree of the alignments between the proto-halo spin directions and the principal axes 
of the local tidal shear tensors would not be so high. 



Motivated by this numerical clue, iLee &: PenI (120001 ) suggested the following quadratic 
formula to quantify the expected degree of the alignments between the halo spins and the 
local tidal shears, generalizing the linear tidal torque theory: 



1 + a 



6ij — aTikTi 



ik -'- kj ■ 



(3) 



where (Lj) is the unit spin vector of a halo, (Tjj) is the unit traceless tidal tensor smoothed on 
the halo mass scale, and a is a correlation parameter in the range of [0, 3/5] which measures 
the strength of the correlation between (Lj) and (Tij). It is worth recalling the fact that 
the ensemble average (LiLj) in the left hand side of equation ([3]) is obtained by taking the 
average of LiLj from the sample halos having all different mass, while the tidal shear field in 
the right hand side of equation ([3]) is smoothed on the single mass scale which amounts to 
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the mean halo mass (i.e., the mean value of the mass of the sample halos) ]Lee fc Pen! (120001 ) 
derived equation ([3]) empirically, relating the angular momenta of the halos with different 
mass to the tidal shear field smoothed on the mean halo mass. 

The correlation parameter a in equation ([3]) which characterizes equation ([3]) is intro- 
duced to take into account strong correlations between (Jjj) and (Tjj) and any modification 
in the subsequent evolution after the moment of the turn-around. The two extreme cases of 
a = and a = 3/5 represent no correlation and maximum correlation between (Lj) and (Tj^), 
respectively. It is expected that the true value of a would be between these two extremes. 

The linear tidal torque theory itself cannot make any prediction on the value of a. Thus, 
the true value of this free parameter a has to be determined empirically. So far, several 
attempts have bee n made to measure the value of a e i ther from numerical simulations or 



from observations (ILee &: Penll2000l : iPorciani et al.ll2002l : iNavarro et. al.ll2004j : iTrujillo et. al. 



20061 : iLee fc Erdogdul 120071 ). Although no general consensus has been reached on the true 



value of a, most of the numerical and observational evidences indicated that the true value 
of a deviates from zero. 



2.2. The First Order Approximation to the Halo Spin Correlations 

A difficulty in constraining the value of a lies in the fact that it is hard to measure the 
principal axes of the linear tidal field in real space. An alternative app roach to the tru e 



value of a is to measure the intrinsic spin-spin correlations of dark halos (jPen et al.l 120001 ). 



The local alignments between the halo spin vectors and the tidal shear te nsors will l ead to 



spatial spin alignments between the neighbor halos. Using equation iPen et al.l (120001 ) 
derived the following analytic model for the halo spin-spin correlations: 

r/(r) ^ (I J(x) ■ J(x + r)H - i ^ ^Pl^" 

Here, the constant 1/3 represents the value of (|J(x) ■ J(x + r)p) for the case of no alignment, 
and ^(r; R) is the linear two-point correlation function defined as 

^{r-R) = r A^{k)^^W\kR)d\nk, (5) 

J-oo kr 

where A^(fc) is the dimensionless power spectrum and W{kR) is the top-hat spherical filter 
of scale radius, R, which is related to the halo mass scale M asR= [?>M / {A'Kp)]^/^ with the 
mean density p. 
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Note that equation (jlj) is proportional to the square of the two point density correlation 
^(r), which means that the halo spin-spin correlations would decrease very rapidly. For in- 
stance, on a galactic scale {R ~ lh~^Mpc) , it is expected that the halo spin-spin correlations 
exist only at distances of order of a few Mpcs and effectively vanishes at larger distances. 

Since the tidal field on different scales are cross-correlated, equation ([3]) leads to the 
existence of the cross-correlations of the spin axes of halos on different mass scales. In 
accordance with equation (jl]) the halo spin cross-correlations on different scales and R2) 
are modeled as 

where ic{r] Ri, R2) is the density cross-correlation function defined as 

A\k)——W{kRi)W{kR2)d\nk. (7) 
■00 kr 

It is worth mentioning that the value of a for the cross correlation ric{r) would not be same as 
its value for the auto-correlation T]{r) in equation (jlj) due to the difference in the smoothing 
scale for the tidal field. Note that equation ([6]) is a modified version of equation (jlj), both of 
which are based on equation ([3]). The correlation parameter a in equation (jlj) represents the 
spin-spin correlations caused by the spatial correlations of the tidal fields smoothed on the 
same scale (one mean halo mass scale, R). While the correlation parameter a in equation (6) 
represents the spin-spin correlations caused by the spatial correlations of the tidal fields on 
two different scales (two mean halo mass scales, Ri and i?2)- Since the spatial correlations 
of the two tidal fields smoothed on two different scales, Ri and R2 (say, -Ri < R2), are lower 
than the case that the two tidal fields are smoothed on the same scale of R2 and higher than 
the case that they are smoothed on the same scale of -Ri, it is naturally expected that the 
value of a for the cross-correlation ric{r) should be different from that for the correlation 
ri{r). 



2.3. The Nonlinear Effect and Large Scale Correlations of Halo Spins 



As mentioned in §1, it was lHui fc Zhand (120021 ) who first pointed out that equation (jl]) is 
valid only if the gravitational tidal field is Gaussian. They argued that the non-Gaussianity 
in the tidal field should not be negligible in the nonlinear regime and suggested that the 
halo spin-spin alignments should be better described as a linear scaling with the density 
correlation function on sufficiently large scales. 



Inspired by the work of 



Hui &: ZhangI (120021 ). here we suggest a new formula to model 
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the halo spin-spin correlations: 

where a nonlinear correlation parameter is introduced to measure the strength of the 
nonlinear tidal effect on the halo spin-spin correlations. In this model the linear correlation 
parameter is notated as ai since its value could be different between the two cases that the 
nonlinear effect is ignored as in equation (jl]) and taken into account as in ([8]). 

Equation ([8]) implies that if the nonlinear contribution exists, then the halo spin-spin 
correlations must exist on larger scale than predicted by the linear model (eq. [4]) since 
it scales linearly with the density correlation function in the nonlinear model. That is, 
the nonlinear tidal effect would generate large-scale correlations of the halo spin axes. Of 
course, the values of the nonlinear correlation parameter as well as the linear correlation 
parameter have to be determined empirically. 

The cross-correlation, ric{r) can be also modeled in accordance with ([H]) as 

' ^ 6 eS(0; Ru R2) ^ "^ec(0; Ru R2) ' ^ ' 

It is worth mentioning that the value of Cni for the cross correlation ric{r) would not be 
necessarily lower than its value for the auto-correlation ri{r) in equation (IHl), unlike the case 
of the linear tidal torque model in §2.2. In the nonlinear regime where the non-Gaussianity in 
the density field grows, the occurrence of the halo merging is quite frequent. This occurrence 
of the halo merging plays a role of transferring the orbital angular momentum generated by 
the external tidal field into the spin angular momentum of a merged halo. In other words, the 
orbital angular momentum of a system composed of small individual halos become the spin 
angular momentum of a large halo formed through merging of the small halos. Given that the 
orbital angular mom entum of the system b efore merging is generated by the external tidal 



field on larger scale (jVitvitska et. al.ll2002l ). the transfer of the orbital angular momentum 
into the spin angular momentum creates cross-correlations of the tidal fields between different 
scales, which in turn generates large-scale cross-correlations rjc in the nonlinear regime. 
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3. SIGNALS FROM NUMERICAL SIMULATIONS 

3.1. Numerical Data 

We use the halo catalogs from the Millennium Run Simulation of 10^° dark matter 
particles for a ACDM uni verse with the cosm ological parameters given as Q^n = 0.25, h = 



0.73, cTg = 0.9 and Us = 1 (jSpringel et al.ll2005l ). The simulation was performed in a periodic 



box of linear size of 500h ^Mpc to follow the evolution of all particles from z = 200 to z = 0, 
each of which has a mass of 8.6 x IO^H'^Mq. 

We analyze the halo catalogs at four different redshifts z = 0, 0.5, 1 and 2 which contain 
information on the halo mass, comoving position, comoving peculiar velocity and the spin 
vector. At each redshift, we restrict our attention only to those halos which have more than 
200 particles for better accuracy in the measurement of the halo spin vector. It amounts to 
selecting only halos which are more massive than 17.2 x IO^H'^Mq. 

Basically, we measure the spatial correlations of the spin axes between neighbor halos 
as a function of separation distance and investigate how the correlation strength changes 
with redshift (z), halo mass (M), velocity [v) and specific angular momentum (/, angular 
momentum per mass). Then, we compare the numerical results with the linear and the 
nonlinear analytic models described in §2 to determine the best-fit values of the correlation 
parameters, a, ai, and 6^1- 



3.2. Dependence on the Redshift 

Figure [T] plots the halo spin-spin correlation function r]{r) at z = 0,0.5,1 and 2 in the 
top- left, top-right, bottom- left, and bottom- right panel, respectively. In each panel, the solid 
dots correspond to the numerical results from the Millennium data while the solid and the 
dashed lines represent the nonlinear model (eq.[8]) and the linear model (eq. [1]), respectively. 
The case of no alignment is also shown as a dotted line in each panel for comparison. 

The errors of the numerical results are calculated as one standard deviation for the case 
of no alignment. Basically, it is computed by the formula of 2/\/45n/^ where is the number 



of halos belonging to each bin (ILee fc Peru l200lh . Note that the sizes of the vertical errors 



in Figure [H depend on the value of Uh which in turn depends on the bin size of the spatial 
separation r. Here we choose the bin size of r as 2/i~^Mpc for all plots. We have tested 



It is now available at 



littp://www. mpa-garching.mpg.de/millennium 
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whether or not the final results depend on the bin size of r and confirmed the robustness of 
the final results. 

For the analytic models, the correlation parameters are determined through fitting to 
the numerical results by means of the x^-minimization. For the eval uation of the analytic 



models, we use the formula for the linear power spectrum given by iBardeen et al.l (119861 ) 
using the same values of the cosmological parameters that are used for the Millennium Run 
simulations and the shape parameter F = Qmh (private communication with V. Springel). 

Table [1] lists the number of halos (Nh), the mean mass (M), the best-fit values of a, a^i 
and ^ni at 2; = — 2. Here, the mean mass M for each bin is used to find the Lagrangian 
smoothing scale, R which has to be implemented into the analytic models (eqsJH -|8]). As 
shown in Figure (H we detect clear signals of the halo spin-spin alignments within a few 
/i~^Mpc distances at all four redshifts. Note, however, a conspicuous difference between the 
results at low redshifts (2; = and 0.5) and at higher redshifts {z = 1 and 2). At 2; = 1 
and 2, the correlation signals rapidly decrease with distance and disappear at distances 
around 10/i~^Mpc. Thus, for the low redshifts, the linear model (dashed line) for itself fits 
the numerical data pretty well. Whereas, at 2 = and 0.5 the correlation signals decrease 
rather slowly, still existent to a nennegligible level at large separations beyond lO/i^^Mpc. 
Thus, at low redshifts the nonlinear model (solid line) fits the numerical data much better 
than the linear model, predicting the existence of the large-scale correlations. 

This phenomena can be quantified in terms of the bestfit-values of the correlation pa- 
rameters listed in Tabled! At 2, the best-fit value of ^ni is effectively zero within one standard 
deviation (cr) while the best-fit value of a\ is higher than lOa (not marginalized level). In 
contrast, at 2; = the value of ^ni deviates from zero as significantly as 12(T while the value 
of a\ decreases significantly to a 3cr level. It is, however, worthing mentioning here that 
the standard deviation, a, here is not marginalized. Figure [2] plots the 68%, 95% and 99% 
contours for a\ and ^ni at ;z = — 2, showing the degeneracy in putting the constraints on 
the two parameters. As can be seen, at ;s = and 0.5, the values of e^i deviate from zero at 
higher than 99% level. 

Our numerical detection implies that the linear model for the halo spin-spin correlations 
(eq.[l]) works well at 2 > 1 but the nonlinear effect indeed dominates at low redshifts 
(2; < 1), producing large-scale correlations of the halo spins, which can be well described by 
our nonlinear model (eq. [8]). 
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3.3. Dependence on the Specific Angular Momentum Magnitude 

We define tlie magnitude of tlie rescaled specific angular momentum (angular momentum 
per unit mass) as / = L/M* where L and M^, are the halo angular momentum magnitude and 
the halo mass in unit of lQ^^h~^MQ, respectively. To determine quantitatively how the halo 
spin-spin correlation changes with Z, we classify the halos at z = into four logarithmic bins of 
/ and measure the correlations separately for each bin. Figure [3] plots 77 (r) at four logarithmic 
bins of / in the top-left, top- right, bottom- left, and bottom-right panels, respectively. Table 
[2] also lists the number of halos, the mean mass, the best-fit values of a, a\ and ^ni for each 
bin, and Figure H] plots the 68%, 95% and 99% contours for a\ andeni at the four bins. 

As shown in Figure [3l in the bins of low-Z halos (top two panels), the correlations are 
rather weak and rapidly diminish to zero at separations greater than a few /;,~^Mpc. Thus, 
for this case the linear model with the single parameter a fits the data pretty well. In other 
words, for the spin-spin correlations of low-/ halos, the values of are effectively zero as 
shown in Table [21 While in the bins of high-/ halos, the correlations are stronger, still existent 
even at separations as large as 10/i~^Mpc. Therefore for the spin-spin correlations of high-/ 
halos, the nonlinear model with non-zero value of ^ni fits the data much better than the 
linear model. Note that for the case of highest-/ halos, the value of ^ni is approximately 
0.005, which is five times greater than the average value of 0.001 at z = 0. 

This phenomena that the halos with higher specific angular momentum tend to have 
stronger non-linear (and linear) spin-spin correlations may be explained as follows. If a halo 
is located in a region where the effect of the tidal field is very strong, then the halo will 
possess high specific angular momentum. Furthermore, its spin direction will be strongly 
correlated with that of its neighbors since the strong effect of the tidal field tends to diminish 
any randomization of the spin orientations that can occur in the subsequent stages. Therefore, 
those halos which have higher specific angular momentum are likely to show strong spin-spin 
correlations under the strong effect of the tidal field. We have also explored whether this 
phenomena occur at higher redshifts z = 0.5,1 and 2 and found that this is still the case. 

3.4. Dependence on the Mass Scale 

The halos at 2; = are classified into four bins according to the value of the rescaled 
mass and their spin-spin alignments are measured separately for each bin. Figure [5] plots 
r]{r), Table [3] lists the number of halos, the mean mass, the best-fit values of a, ai and e^\, 
and Figure [2] plots the 68%, 95% and 99% contours for ai and ^ni for the four bins. 

As shown in Figure O in the third bin which includes the galactic halos with mass of 
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order of (1 — 2) x 10^^/i~^Mq (bottom-left panel), the correlation signal is strongest, existent 
at separations as large as lO/i'^Mpc. For this bin, the nonlinear model with nonzero value 
of ^ni gives a better fit the numerical results than the linear with model single parameter a. 
The best-fit value of ^ni deviates from zero by a factor of Sex as shown in Table [31 Whereas in 
the first and the fourth bins which include dwarf halos and galaxy-group halos, respectively, 
the signals are rather weak and the correlations disappear at large distances. For these 
bins, the linear model for itself fits the numerical data quite well and the values of ^ni are 
effectively zero. This result implies that the strength of the nonlinear effect depends on mass 
scales, being strongest on the galactic mass scale. Note, however, here that the value of the 
characteristic mass scale changes with z. Thus, the results given in §3.2 and §3.4 are related 
to each other. 



3.5. Scale Bias and Cross-Correlations 

As mentioned in §2.3, if the nonlinear effect dominates on the halo spin-spin correlations, 
then it is likely to produce stronger cross-correlations between different mass scales (eq. [9]). 
We measure the halo spin cross-correlations between the halos of different mass bins that 
are constructed in §3.2 (the first bin corresponds to the lowest mass scales, see Table [3]), and 
find the best-fit values of a, a„i and ^ni by fitting the numerical data to the analytic model 
(eq.il). 

Figure [7] plots the numerical results of rjc for the six different cases: 77ci2, tlciZi Vcm, 
Vc23, Vc24:,Vc34:, where r]cij represents the cross-correlations between the i-th and the j-th 
mass bins. As can be seen, the nonlinear model basically gives better fit to the numerical 
result for all cases, indicating that the nonlinear effect increases indeed increase spin cross- 
correlations between different mass scales, as presumed in §2.3. Figure [9] plots the 68%, 95% 
and 99% contours for a\ and ^ni for these two cases of r/is and 7723 where the strong signals of 
the cross-correlations are detected. Obviously, the values of e^i for each case deviates from 
zero at 95% confidence level. 



3.6. Dependence on the Velocity Magnitude 

We also investigate the dependence of the halo spin-spin correlations on the magnitude 
of halo peculiar velocity v in unit of km/s. The correlations rj{r) are plotted in Fig{TO]at four 
logarithmic bins of v, and the 68%, 95% and 99% contours for the correlation parameters of 
the nonlinear models are shown in FigJTTl 
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As shown in Figure [TOl the dependence of the halo spin-spin correlations on the velocity 
magnitude turns out not to be so strong. The significant nonlinear effect on the halo-halo 
correlation is found only in the first bin of low-f halos (top- left panel), For this bin, the 
linear model fails in fitting the numerical data while the nonlinear model gives a better fit, 
predicting the large-scale correlations. As shown in the top-left panel of FigJTT| the value of 
£ni for the case of lowest-f halos deviates from zero at 95% confidence level. 



4. SUMMARY AND DISCUSSION 

We summarize our results in the following: 

• To account for the possible nonlinear tidal effect on the intrinsic halo alignments, 
we suggest a new model characterized by two free parameters, oi and e^i- The two 
parameters ai and represents the strength of the linear and the nonlinear tidal 
effect on the orientations of the halo spins, respectively. If ^ni = 0, it corresponds to 
the case that the nonlinear tidal effect is negligible and the halo spin-spin correlations 
is a quadratic scaling of the density correlation function ^(r), existent on small scales. 
While if £ni deviates from zero, the nonlinear tidal effect is nennegligible and the halo 
spin-spin correlations is a linear scaling of C,{r), present on large distance scales. 

• Using the halo catalogs from numerical simulations at z = 0,1,0.5 and 2, we have 
measured the halo spin-spin correlations and determined the best-fit values of oi and 
£ni- It is found that the values of ai increases with z while e-^i decreases with z. At 
z = 1 and 2, the values of are found to be negligibly small, while at z = 0.5 and 
z = Erd has an order of 10~^ higher than 5a. Especially at z = 0, ai is effectively zero 
and Eni deviates from zero by a factor more than lOcr, which implies that at present 
epoch the nonlinear tidal effect is dominant inducing halo spin correlations on scales 
as large as lO/i^^Mpc. 

• We have investigated the dependence of a\ and ^ni on the halo intrinsic properties such 
as mass (M), specific angular momentum magnitude (Z), and velocity magnitude (v) at 
z = 0. It is found that the value of ^ni depends most sensitively on the value of I. For 
those halos with highest /, the value of ^ni reaches as high as 5 x 10^^ and the value of 
ai is effectively zero. Regarding the dependence on M and v, it is found that ei is less 
sensitive. However, for those halos with M ~ (1 — 2) x 10^^/i~^Mq logf < 2.5km/s, 
the value of ^ni tend to be higher. 

• We have measured the cross-correlations of the spin axes between halos on different 
mass scales. It is found that on two different galactic mass scales (Mi ~ IO^^H^^Mq 
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and Ms ~ 1O^^/;,~^M0, the cross-correlations exist at scales as large as 10/i~^Mpc, due 
to the growth of the non-Gaussianity increases cross-correlations of the tidal field on 
different scales. 



Our results have shown that the galaxy intrinsic alignments can exist not only for the 
case that the separation distance between the two galaxies is small but also for the case that 
any galaxy pair separated by distances larger than lO/i'^Mpc can still have aligned axes. It 
also implies pr eferential E-mode conta mination of the galaxy intrinsic alignments in weak 
lensing signals (ICrittenden et al.ll2002l ). 



In addition, we have also found quantitatively that the galaxy intrinsic alignments could 
be cross-correlated with the larger-scale tidal field. These results imply that the intrinsic 
galaxy correlations are indeed correlated with the surrounding density field on large scales, 
which will in turn lead to correlation between the gravitational lensing shears and the intrinsic 
galaxy alignments, (often called, GI correlations). 



Recently, iHirata et al.l (120071 ) claimed that a signal of the GI correlations has bee n 
detected from the large low-redshift galaxy surveys (see also, iMandelbaum et al.l l2006l ). 
According to their claim, the GI correlations exist to a nennegligible level for the bright 
Luminous Red Galaxies (LRGs) and it is zero for the faint blue galaxies, mentioning an 
uncertainty in the estimation of the GI correlation for the faint blue galaxies. According 
to our numerical detection, the nonlinear tidal effect which causes the large-scale intrinsic 
alignments is strongest for the halos with high specific angular momentum magnitude and 
l ow peculiar velocity , which usually correspond to the galaxi es with low surfac e brightness 
( iJimenez et allll998l . and references therein) and blue color (jPavis et al.l 119971 ). Thus, our 
results indicate that the GI correlati on amplitude f or the faint blue galaxies, may exceed the 
estimation upper limit, as noted by iHirata et al.l (120071 ). It will be of interest to compare 
the cross-correlations between the spin axes of the blue galaxies and the shapes of the red 
galaxies for the comparison of their correlation scales. 



The Millennium Run simulation used in this paper was carried out by the Virgo Super- 
computing Consortium at the Computing Centre of the Max-Planck Society in Garching. 
We thank V. Springel and G. Lemson for useful comments. We also thank an anonymous 
referee for helpful suggestions. J.L. is very grateful to the warm hospitality of L. Kofman 
and the Canadian Institute for Theoretical Astrophysics (CITA) where this work was initi- 
ated and performed. J.L. acknowledges the financial support from the Korea Science and 
Engineering Foundation (KOSEF) grant funded by the Korean Government (MOST, NO. 
ROl-2007-000-10246-0). 
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Table 1. The number of halos, the mean mass, the best-fit values of the hnear and the 
nonhnear correlations parameters for the intrinsic halo spin correlations at four redshifts 

from the Millennium Run data. 



z 




(IO^o/i-^Mq) 


a X 10^ 


ax X 10^ 


£nl X 10^ 





1846776 


169.7 


11.3 ±0.6 


3.9 ± 1.6 


1.2 ±0.1 


0.5 


1942263 


133.3 


9.8 ±0.7 


6.1 ± 1.1 


0.7 ±0.1 


1 


1928690 


106.7 


8.6 ±0.8 


7.9 ± 0.9 


0.1 ±0.1 


2 


1558504 


69.6 


10.6 ±0.9 


10.4 ±0.9 


0.0 ±0.2 
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Table 2. Numerical results for the intrinsic halo spin correlations at four bins of the 

specific angular momentum magnitude /. 





log(0 




(lO^o/i-iM©) 


a X 102 


ai X 102 


£nl X 10^ 


(■ 


-4.0, -2.0) 


350541 


607.1 


9.4 ±2.4 


9.4 ±2.4 


0.0 ±0.4 


(■ 


-2.0,-1.7) 


604356 


102.7 


10.1 ±2.2 


6.6 ±3.4 


0.7 ±0.4 


(- 


-1.7,-1.5) 


475842 


51.0 


16.8 ±2.1 


5.4 ±6.3 


2.7 ±0.6 


( 


-1.5,0.0) 


416032 


32.5 


20.2 ±2.2 


0.0 ±9.7 


4.9 ±0.8 
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Table 3. Numerical results for the intrinsic halo spin correlations at four different mass 

ranges. 









a X 10^ 


ai X 10^ 


£nl X 10^ 


(lO^o/i-^Mo) 




(lO^o/i-^M©) 








(17.2,38.0) 


925011 


25.0 


14.1 ±1.6 


8.5 ±2.6 


1.3 ±0.4 


(38.0,85.0) 


469647 


55.7 


14.1 ±2.3 


10.6 ±3.2 


1.0 ±0.6 


(85.0,190.0) 


233724 


124.6 


20.8 ±2.5 


15.0 ±3.5 


2.4 ±0.9 


(190.0, 10^) 


218248 


1022.5 


15.5 ± 1.9 


13.0 ±2.2 


0.8 ±0.6 
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i 



1 10 1 10 

r [h"iMpc] r [h-^Mpc] 

Fig. 1. — Intrinsic halo spin-spin correlations at four different redshifts: z = 0, 0.5, 1 and 
2 (top-left, top-right, bottom-left, and bottom-right, respectively). In each panel, the solid 
dots represent the numerical results, while the solid and dashed lines correspond to the 
nonlinear (eq. [1]) and linear (eq. [1]) analytic model, respectively. The errors are calculated 
as the standard deviation for the case of no correlation. The dotted line represents the case 
of no correlation. 




Fig. 2. — 68%, 95% and 99% contour plots for the linear and the nonlinear correlation 
parameters, a\ and £ni at 2; = 0, 0.5, 1 and 2 (top- left, top-right, bottom- left, and bottom- 
right, respectively). In each panel, the cross mark represents the best- fit value. 
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Fig. 3. — Intrinsic halo spin-spin correlations ai z — for the four different logarithmic bins 
of the specific angular momentum magnitude Z at 2; = 0. 




Fig. 4.— 



68%, 95% and 99% contour plots for ai and e^iat the four logarithmic bins of /. 
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Fig. 5. — The dark halo spin-spin correlations at 2; = for the four different bins of the halo 
mass in unit of 10^°/i-^Mq. 




Fig. 6.— 



68%, 95% and 99% contour plots for ai and at the four bins of halo mass M^. 
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1 10 1 10 

r [h"iMpc] r [h-^Mpc] 

Fig. 7. — Intrinsic halo spin cross correlations between different mass scales at z = 0. In each 
panel, rjcij represents the cross-correlation between the i-th and the j-th mass bin, given in 
Table [31 
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Fig. 8. — Comparison of the spin correlation r]ii{r) between halos in the same mass bins with 
the spin cross-correlation rjajir) between halos in different mass bins. 
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Fig. 9. — 68%, 95% and 99% contour plots for ai and £ni for the spin cross-correlations 
between halos in the 1st mass bin and the 2nd mass bin (top) and between halos in the 2nd 
mass bin and the 3rd mass bin (bottom). 
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Fig. 10. — Intrinsic halo spin correlations at the four logarithmic bins of the velocity mag- 
nitude V in unit of km/s. 
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Fig. 11.— 



68%, 95% and 99% contour plots for ai and £ni at the four logarithmic bins of v. 



